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Abstract: The idea of a variable dark energy has been entertained many times in the 
literature and from many different points of view. Quintessence is just a popular way to 
implement this idea in recent times, but so far with little success. Another possibility 
is to think of the cosmological term, A, as a "running quantity" much in the same way 
as the electromagnetic coupling constant. However, the fact that A is a dimension-four 
parameter implies that it may obey a peculiar renormalization group equation, which 
at low energies could be dominated by "soft decoupling" contributions of the form A ~ 
H 2 Alp stemming from physics near the Planck scale. This value lies in the ballpark of 
the measurements from CMB and high-z supernovae. A "renormalized" FLRW cosmology 
of this kind may reveal itself as a sound, and testable, proposal for a variable A model 
within quantum field theory in curved space time. 



1. Introduction 

The Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological model ]l]] constitutes 
the standard paradigm of present day cosmology. Its 4-curvature is determined from the 
various contributions to its total energy-momentum tensor density, namely in the form of 
matter energy density, radiation pressure and cosmological constant: 

R = 8ttG n (p- 3p + 4A) . (1.1) 

The cosmological constant contribution to the curvature of space-time is represented by 
the A term. The latter enters the original gravitational field equations in the form 

Rfiv - -^9nvR = SttGn Tju> , (1-2) 

'Speaker. 

t Associated with Instituto de Ciencias del Espacio-CSIC. 
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where T^ u is given by T^ u = + g^A^t), T^ u being the ordinary energy-momentum 
tensor associated to isotropic matter and radiation. Here we entertain the possibility that 
A is a function of time 1 . By the Bianchi identities, it follows that A is a constant if and 
only if the ordinary energy-momentum tensor is individually conserved {xj^T^ = 0). In 
particular, A must be a constant if is zero (e.g. during inflation). 

Modeling the expanding universe as a perfect fluid with velocity 4- vector field U^, we 
have 

T^ = -pg^ + {p + p)U»U\ (1.3) 

where p is the isotropic pressure and p is the proper energy density of matter. Clearly the 
modified T^ v defined above takes the same form as ([D|) with p — > p = p + A , p — > p = 
p — A. With this generalized energy-momentum tensor, and in the FLRW metric (k = 
for flat, k = ±1 for spatially curved, universes) 

ds 2 = dt 2 - a 2 (t) — r + r 2 d6 2 + r 2 sin 2 9 d<)) 2 ) , (1.4) 

\ 1 — k r z J 

the gravitational field equations boil down to the Friedmann-Lemaitre equation 



" I 1 ) 2 — ^!/;-A) -• (1.5) 



< ^) =^(p + A)-4 
a J 3 

and the dynamical field equation for the scale factor: 



4tt Air 

-—G N {p + 3p)a = -—G N (p + 3p-2A)a. (1.6) 



The physical CC is the parameter A involved in the equations ( |1.5| ) and ( |1.6[) . As these 
equations are used to fit the cosmological data ||, |3| , the parameter A involved in them is 
to be considered the physical (observed) value of the cosmological constant. A first integral 
of the system ( |1.5[ ) and ( |1.6| ) is given by 

A + p + 3H{p + p) = 0, (1.7) 

which shows that a time-variable CC cosmology allows transfer of energy from matter- 
radiation into vacuum energy, and vice versa. Since the radiation pressure is negligible 
at present, Eq. ( |1.5| ) can be rewritten as the following exact sum rule for the cosmological 
parameters: 

n° M + n° A + n° K = i, (1.8) 



with 



Pc Pc 11 a 



p® M and Ao being the matter density and cosmological constant at our time, and 
_3Hi_ 

8 7T G/V 

1 While the homogeneous and isotropic FLRW cosmologies do not allow spatial gradients of A, they do 
not forbid the possibility that A may be a function of the cosmological time: A = A(t). 



p° c = ~ 10.6 hi GeV/m 3 ~ (3.0 y% x 10^ 3 evY (1.10) 
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is the present value of the critical density (equivalent to a few protons per cubic meter). 
Here the dimensionless number ho ~ 0.7 ± 0.1 sets the typical range for today's value of 
Hubble's constant Hq = (d/a) = 100 (Km/ sec Mpc) Hq. 

What is the nature of the total energy density of the universe? On the one hand the 
baryonic contribution to the total matter content pu (a scanty £Ib ~ 5% of the critical 
density, derived from nucleosynthesis calculations) is far smaller than the total amount of 
matter detected by dynamical means |l], [|, ||, namely ~ 30% of the critical density. 
Therefore the bulk of the matter content of the universe is dark matter, and must be in the 
form of an unknown kind of cold (non-relativistic and non-baryonic) invisible component. 
Significant amounts of hot (relativistic) dark matter are excluded because it would not 
fit with the models of structure formation (particularly at small scales). So the radiation 
part at present boils down to an insignificant (few per mil) fraction of neutrinos (yet 
comparable to the total amount of luminous matter that we see!) plus an even more 
negligible contribution of very soft photons (the relic CMB radiation) entering at the level 
of one ten-thousandth, at most, of the critical density today ||. 

On the other hand, the astrophysical measurements tracing the rate of expansion of the 
universe with high-z Type la supernovae ||, [| indicate that ~ 70% of the critical energy 
density of the universe is cosmological constant (CC) or a dark energy candidate with a 
similar dynamical impact on the evolution of the expansion of the universe. Specifically, 
the CC value found from Type la supernovae at high z is: 

A = p l ~ 6 hi x 10~ 47 GeV 4 . (1.11) 

Independent from these supernovae measurements, the CMB anisotropies Q, including the 
recent data from the WMAP satellite, lead to = 1.02 ± 0.02 B. As a first observation, 
it is obvious that this result leaves little room for our universe to be spatially curved 
(Qk < 2%), and indeed it suggests that our universe is spatially flat, as expected from 
inflation H . As a second observation, when combining this result with the dynamically 
determined value (from clusters of galaxies) of the matter density (viz. f2j{ ~ 30%), the 
bookkeeping in Eq. (|1.8|) leads us to an outstanding conclusion: the rest of the present 
energy budget (a large gap of order 70% of the critical density) must be encoded in the 
parameter Q^. Hence the CMB measurements and the high-z supernovae data lead to 
a conclusion in concordance with the value of A, if one accepts the data on from 



clusters. As mentioned, the value ( |1.11 ) corresponds to the parameter entering the classical 
cosmological equations ( |1.5| ) and ( |1.6| ) . Therefore the situation seems to be consistent from 
the experimental point of view 2 

What about the theoretical situation? In Quantum Field Theory (QFT) we have long 
expected that the vacuum fluctuations should induce a non- vanishing value for A ||] , and 
the question is whether in realistic QFT's we have a prediction for Ao in the ballpark of 
the measured value ( |l.ll| ). Sadly, the answer is no. For, in the context of the Standard 
Model (SM) of electroweak interactions, this measured CC should be the sum of the original 
vacuum CC in Einstein's equations, A vac , and the induced contribution from the vacuum 



2 See, however, Ref. M for a more critical point of view on the present "Concordance Model" 
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energy of the quantum fields: 

A = A mc + A ind . (1.12) 

What is the expected value for Aj n d in the SM? Let $ be the quantum Higgs field operator. 
In the ground (vacuum) state of the SM, the vacuum expectation value (VEV) of <J> + $ 
will be denoted < <& +< I> >= \<ft 2 , where (ft is the mean field associated to the field operator 
<!>. The classical potential associated to the mean field reads 

V c i = -\ml^ + ^. (1.13) 

Shifting the original field (ft — > H° + v such that the physical scalar field H° has zero VEV 
one obtains the physical mass of the Higgs boson: Mh = y/2 ms- Minimization of the 
potential (1.13) yields the spontaneous symmetry-breaking relation: 



/2m5 M 2 
(ft=\—*=v and A = — f. (1.14) 
V A v z 

The VEV < $ >= v/V2 gives masses to fermions and weak gauge bosons through 

m fi = hi-j=, M 2 v = ±g 2 v 2 , M 2 z = ^(g 2 +g' 2 )v 2 , (1.15) 

where hi are the corresponding Yukawa couplings, and g and g' are the SU(2)l and 
U(1)y gauge couplings. The VEV can be written entirely in terms of the Fermi scale 
Mp = G F 1/2 ~ 293 GeV as follows: v = 2~ 1 /' i M F ~ 246 GeV . From (|LT^) one obtains 
the following value for the potential, at the tree- level, that goes over to the induced CC: 

m \ loo 

Knd=<V cl >=-^ = --M 2 I v 2 . (1.16) 

From the current LEP 200 numerical bound on the Higgs boson mass, Mjj > 114.1 GeV 
||, one finds \A ind \ > 1.0 x 10 8 GeV A . Clearly, |Aj n ^| is 55 orders of magnitude larger than 
the observed CC value ( 1.1 1| ) . Moreover, the Higgs potential gets renormalized at higher 



order in perturbation theory, and therefore it is the value of the effective Higgs potential 
V e ff what matters at the quantum level. These quantum corrections by themselves are 



already much larger than (1.11). Finally, we also note that in general the induced term 
may also get contributions from strong interactions, the so-called quark and gluon vacuum 
condensates. These are also huge as compared to ( |1 . 1 1| ) , but are much smaller than the 



electroweak contribution ( 1.16j ). 



Such discrepancy, the so-called "old" cosmological constant problem (CCP) [|10], 11 
manifests itself in the necessity of enforcing an unnaturally exact fine tuning of the original 
cosmological term A vac in the vacuum action that has to cancel the induced counterpart 
Ajnd within a precision (in the SM) of one part in 10 55 . This big conundrum has triggered 
many theoretical proposals. On the first place there is the longstanding idea of identifying 



the dark energy component with a dynamical scalar field [12, 13]. More recently this 



approach took the popular form of a "quintessence" field slow-rolling down its potential [14] 
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and variations thereof [15]. The main advantage of the quintessence models is that they 
could explain the possibility of an evolving vacuum energy. This may become important in 
case such evolution will be someday detected in the observations, although it is not obvious 
that the models can be easily discriminated [16] . Furthermore, a plethora of suggestions 



came along with string theory developments JL7| and anthropic scenarios |l{| . Other recent 
ideas have been put forward, like the intriguing proposal of non-point-like gravitons at 
sub-millimeter distances or the suggestion of having multiple degenerate vacua [ 20 1 . 



There are however other (less exotic) possibilities, which should be taken into account. 
In a series of recent papers f|l], the idea has been put forward that already in standard 
QFT it should not make much sense to think of the CC as a constant, even if taken as a 
parameter at the classical level, because the Renormalization Group (RG) effects may shift 
away the prescribed value, in particular if the latter is assumed to be zero. Thus, in the RG 
approach one takes a point of view very different from the quintessence proposal, as we deal 
all the time with a "true" cosmological term. It is however a variable one, and therefore 
a time-evolving, or redshift dependent: A = A(z). Although we do not have a QFT of 
gravity where the running of the gravitational and cosmological constants could ultimately 
be substantiated, a semiclassical description within the well established formalism of QFT 
in curved space-time (see e.g. [^3], should be a good starting point. Then, by looking 



at the CCP from the RG point of view [25], the CC becomes a scaling parameter whose 



value should be sensitive to the entire energy history of the universe - in a manner not 
essentially different to, say, the electromagnetic coupling constant. 

The canonical form of renormalization group equation (RGE) for the A term at high 
energy is well known - see e.g. However, at low energy decoupling effects of the 

massive particles may change significantly the structure of this RGE, with important phe- 
nomenological consequences. This idea has been elaborated recently by several authors 



from various interesting points of view [£1], 2t , [27]] . It is not easy to achieve a RG model 
where the CC runs smoothly without fine tuning at the present epoch. In Ref. [28, [2^, 30] 
a successful attempt in this direction has been made, which is based on the possible exis- 
tence of physics near the Planck scale. In the following we review the main features and 
implications of this "RG-cosmology" . 

2. Running cosmological constant 

It is well known that the effective action of vacuum takes the form pi, M 



Svac = J d A x^f^g | aiR^p + a 2 i?^ + a 3 R 2 + a±UR - ~^~q — R ~ A «ic} , (2.1) 

and at low energies the only relevant part for phenomenological considerations is the 
Hilbert-Einstein piece 

She = ~ J d A x^g~ ( R + Kac \ ■ (2.2) 

While the phenomenological impact of the higher derivative terms in ( |2.lD is negligible at 
low energy (namely, at current times), the presence of the parameter A vac is as indispensable 
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as any one of these higher derivative terms to achieve a renormalizable QFT in curved space- 
time 3 . Admittedly the vacuum CC itself, A vac , is not the physical (observable) value of 
the cosmological constant. However, its presence is essential to make possible the matching 
of the theoretical value of A in Eq. ( |1.12| ) with the physical measurement Ql.lip . Here we 
do not address the issue of why the two terms (induced and vacuum terms) cancel with 
the desired precision (the "old" CCP). Rather, we take it for granted, and then study the 
possibility that the cosmological constant evolves according to a RGE whose flow crosses 
the present instant of time with the value ( 1.11| ), and we then investigate whether such 



a RGE, when extrapolated back to the past, gives testable and consistent predictions. In 
particular, we verify whether the predicted value of A at the nucleosynthesis time does not 
disturb the basic features of this momentous regime in the history of the universe. Once 
this is guaranteed, we wish to check whether the predictions of the RGE for nearby past 
times depart significantly from the standard FLRW cosmology, and therefore whether we 
can find out signs that can effectively discriminate among the two. Finally, the same RGE 
can give us a prediction for the future or destiny of our universe. 

Of course the first issue is to elucidate the structure of the RGE for the A term at any 
given energy regime. A first clue is to recognize that this renormalization group equation 
is well- know at high energy [24, El]]. For particles i = 1,2,3,... of masses rrii and spins Jj 



one finds |21]: 



where 

Ai = {-l) 2Ji (Ji + 1/2) nji N c , (2.4) 

with 71/0,1,1/2} = (1,1,2) and iV c = 1,3 for uncolored and colored particles respectively. 
However, this equation applies only if the masses m, do correspond to active degrees of 
freedom (d.o.f.), namely if the RG scale ^ satisfies ji 3> m; for all % = 1,2,3,.... Any 
mass that does not satisfy this condition must be dropped from that RGE. At the same 
time this raises two issues: i) One is the meaning of the RG scale [i in cosmology, and the 
other is: ii) what happens with the decoupled contributions, i.e. those heavy masses that 
satisfy M% > fi for a given value of the RG scale. This question has been addressed in 



Ref. 1 21, 26] where it was first recognized that the contribution from the heavy degrees of 
freedom could be relevant, especially in the context of the cosmological constant problem, 
where every "small" contribution (for the standards of high energy physics) can be very 
big (for the standards of cosmology). This reflection suggests that Eq. ( |2.3[ ) can only be 
a "sharp cutoff" approximation to the complete RGE for the CC. In other words, it is 
not enough to (abruptly) couple or decouple the various degrees of freedom following a 
^-function procedure like in the renormalization group equations that we are used to in 



the Minimal Subtraction ( MS )-scheme [plf . Here the smoothing of the various contribu- 
tions could play a role. In fact, this smoothing is perfectly possible if one adopts a more 



3 From this sole observation it follows that quintessence models without a A term cannot be renormaliz- 
able theories in curved space-time. 
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physical renormalization scheme, like the Momentum Subtraction Scheme [^]. It can be of 
course much more cumbersome, but it provides a well-defined prescription to "connect" and 
"disconnect" the various d.o.f. automatically, i.e. not just by hand. Then, from general 
arguments of effective field theory, in combination with the Appelquist-Carazzone (AC) 
decoupling theorem |33| ) , we expect that the generalization of the RGE for any energy 
regime is the following: 

' \ i 3 3 3 3 J 

where the sums are taken over all massive fields; A, B,C, D, ... are constant coefficients, 
and \x is the energy scale associated to the RG running. The r.h.s. of (^5|) defines the 
/?A-function for A, which is a function of the masses and in general also of the ratios of the 
RG scale and the masses, /i/Mj. 

Concerning the interpretation of the RG scale ^ in cosmology, within this, more phys- 
ical, RGE we adopt the proposal of Ref. [^] , namely we identify fi with the value of the 
Hubble parameter at the corresponding epoch: 

H~H(t). (2.6) 

From Eq. (|l.l| ) and ( |l.5| ) it is clear, within order of magnitude, that in the FLRW cos- 
mologies this is a particularization of identifying fi with i? 1 / 2 , the latter being a perfectly 
sensible (and invariant) way to define the RG scale in a general cosmological framework 4 . 
On physical grounds, what this implies is to identify /i in cosmology with the typical energy- 
momentum of the cosmological gravitons for a given metric, in this case the FLRW one 
Ql.4| ). With this Ansatz, it is clear that only particles whose masses are below the value 
of the Hubble parameter (rrii < H) can be active d.o.f. at the given epoch. Notice that 
this makes extremely difficult to find out active particles, from the point of view of the 
RG, in any cosmological epoch! For example, the present value of the Hubble parameter is 
H ~ 1(T 33 eV. The latter is 30 orders of magnitude smaller than the mass of the lightest 
neutrino, 41 orders of magnitude smaller than the QCD scale and 61 orders of magnitude 
smaller than the Planck scale. Obviously, all massive particles decouple the same way! 
Nevertheless, this effective decoupling of all the d.o.f. corresponding to the SM particles 
(in fact of all particles well below the Planck scale!) does not necessary mean that there 



is no running of the CC. The general RGE (|2.5D contains indeed the clue to the possible 



solution of this problem. Namely, A could effectively run mainly due to the heaviest d.o.f., 
represented by the remaining terms on the r.h.s. of ( |2.5D other than those in ( |2.3| ). Look- 
ing at the chain of decoupled terms, we immediately discover the dominant ones in the 
/3A-function, 



4 Scale (2.6) has also been used successfully in other frameworks, e.g. in |34j to describe the decoupling 



of massive particles in anomaly-induced inflation -or modified Starobinsky model [ 
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which we call the "soft decoupling" contributions to the RGE of the CC. These are rather 
peculiar, as they grow quadratically with the masses of the various particles. Notice that 
the existence of these soft decoupling contributions stems from the dimension-four nature 
of the CC. There is no other parameter either in the SM or in the GUT models with such 
property. 

The structure of the RGE (pT5|) is dictated by the the AC-decoupling theorem pi, pi] 



and general covariance. Indeed, dimensional analysis is not enough to explain the most 
general structure of (5\. The fact that only even powers of are involved stems from 
the covariance of the effective action under the identification (|2.6|). Odd powers of \i = 
H cannot appear after integrating out the higher derivative terms as they must appear 
bilinearly in the contractions with the metric tensor. In particular, covariance forbids 
the terms of first order in H. As a result the expansion must start at the i/ 2 -order. In 
short, when applying the AC theorem in its very standard form to the computation of 
the /^-function of the cosmological constant, the decoupling of the heavy masses does still 
introduce inverse power suppression, but since the /?A-function itself is proportional to the 
fourth power of these masses it eventually entails a decoupling law 1/M 2n_4 , and so the 
n = 1 and n = 2 terms do not decouple in the ordinary sense whereas the n = 3 and above 
terms do. The upshot is that the n = 1 (soft decoupling terms) dominate the /^-function. 

What about the phenomenological significance of this RG framework for the cosmolog- 
ical constant? The remarkable thing about the soft decoupling terms is that if we assume 
that there is physics near the Planck scale, characterized by a set of (fermion and boson) 
fields with masses near the Planck scale (Mj < Mp), then we can realize the numerical 
"coincidence" 

cH$M 2 /(4ir) 2 ~ c (1.5 KT 42 GeV x 1.2 10 19 GeV) 2 / (4tt) 2 ~ 1(T 47 GeV 4 ~ A (2.8) 



for some c = 0(1 — 10). Put another way: using the RG scale Ansatz (|2.6| ), the contribution 
from (near-) Planckian size masses to the soft decoupling terms is just of the order of the 



present day value of the cosmological constant - Cf. Eq. (1.11). This means that, thanks 
to the soft decoupling terms fl2.7| ), there is the possibility to obtain an effective running of 
the CC which changes its value by an amount that is of the order of the CC value itself. 
Hence a smooth running of the CC is feasible without requiring any fine tuning of the 
parameters. 

Let us introduce the following mass parameter: 



M = 



I>^ 2 



(2.9) 



The mass Mj of each superheavy particle may be smaller than Mp and the equality, or even 
the effective value M > Mp, can be achieved due to the multiplicities of these particles. 
At the end of the day we expect from fl2.7| ) that the RGE of the CC is dominated by the 
term 

' /A ' aH 2 M 2 + ... = ^H 2 M 2 +... (2.10) 



dhxH (4vr) 2 " 4vr 
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where the signature a = ±1 of the /?A-function depends on whether the fermions (a = — 1) 
or bosons (a = +1) dominate at the highest energies. In the previous RGE we have defined 
the fundamental (dimensionless) parameter 



a M 2 



12 7T M 2 



(2.11) 



which, as we shall see, traces out the presence of the RG effects on the modified cosmology 
as well as its departure from the standard FLRW framework. 

A few reflections are now in order. The physical interpretation of the renormalization 
group in curved space-time requires the formalism beyond the limits of the well-known 
standard techniques. These techniques are essentially based on the minimal subtraction 
MS scheme of renormalization [31, p4|], which does not admit to observe the decoupling. 



The physical interpretation of the renormalization group in the higher derivative sector 
of the vacuum action ( |2.1[ ) can be achieved through the calculation of the polarization 
operator of gravitons arising from the matter loops in linearized gravity [j3"6[| . In this case 
one can perform the calculations in the physical mass-dependent renormalization scheme 
(specifically, in the momentum subtraction scheme, in which \i is traded for an Euclidean 
momentum p), obtain explicit expressions for the physical beta- functions and observe the 
decoupling of massive particles at low energies. In the mass-dependent scheme one has 
direct control of the functional dependence of the /3-functions on the masses of all the 
fields, i.e. one knows explicitly the functions fikiyn^p 2 /Mj) for the higher derivative 
terms in (p.l|). Quite reassuring is the fact that these momentum subtraction (3- 
functions boil down, in the UV limit, to the corresponding /^-functions in the MS scheme 
[36]. At present there is no method of calculations on the non-flat background compatible 



with the physical renormalization scheme. In this situation our phenomenological approach 
looks rather justified. The covariance of the effective action forbids the first order in H 
corrections. Then, recalling that the UV contribution to j3\ from a particle of mass m is 
Pa ~ w 4 , the low-energy contributions to (3\ must be supressed by, at least, the factor 

n/m 2 ~ H 2 /m 2 , (2.12) 



and hence the overall low-energy j3\ acquires the form H 2 m 2 , Eq. (|2.10|) . The same 
form of decoupling can be expected for the induced counterpart of the CC. The 

renormalization group equations for Aj n d and A vac in ( |1.12j ) are independent, and the 
relation emerges only at the moment when we choose the initial point of the RG trajectory 



for the vacuum counterpart [21, [22f| . In fact, the Feynman diagrams corresponding to the 
Ai n d are, at the cosmic scale, very similar to the vacuum ones, because the Higgs field 
<f> is in the vacuum state, < >^ 0, and the non-zero contribution of these diagrams is 
exclusively due to the momenta coming from the graviton external lines. Finally, there is 
no real need to distinguish induced and vacuum CC's in the present context. In both cases 



the absence of the n = order (non-supressed) contributions in (2.5) is required by the 
apparent correctness of the Einstein equations and the smallness of the observable CC. 
Let us finally remark on the running in other sectors. In the basic equation of energy 



conservation (1.7) one could have enforced the conservation of matter by itself, but if one 
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still wishes an evolving CC this could only be posible at the expense of a time-varying 
Newton's constant, which could give rise to various interesting effects |^7|, [37j. However, in 
our framework this is not granted. For the inverse Newton constant 1/G the running is 
irrelevant because 1/G ~ Mp is very large and the effect of the running is relatively small 
as it was demonstrated in plj . On the other hand, the relevance of the higher derivative 
terms at low energies is supposed to be negligible, as it was recently shown in [38] for the 
perturbations of the conformal factor. In contrast, the soft decoupling for the CC does 
matter because the CC is very small and any running, even the very small one that we 
get, is sufficient to produce some measurable effect at large redshifts. For precisely this 
reason we have to use the Friedmann equation, in combination with the full conservation 
law including the variable cosmological term, and no other terms. 



3. A "renormalized" FLRW cosmology with running A 

After introducing and discussing our theoretical framework, let us now solve explicitly for 
the cosmological evolution equations for the matter and vacuum energy densities ensuing 



from our semiclassical FLRW model [ 28 , 30 1 . Let us first of all assume that we are in the 
matter era, then p = 0. We can now solve for the simultaneous system formed by equations 
( |1.5| ), fll.7| ) and our basic RGE ( 2.1C| ). It is convenient to eliminate the time variable in 
favor of the redshift z, using ao/a = 1 + z. The result is obtained after straightforward 
calculation. For the matter density we find 



p{z;u) 



P°M + 



1 Q~Pc 

1 — 6v 



(1 + *) 



3(1-!/) 



1 - 3z/ 



p° c (l + z) 



where we have introduced the parameter 



-2 v£fl K . 



(3.1) 



(3.2) 



It is nonzero only for spatially curved universes, and it depends on the fundamental pa- 



rameter v, Eq. ( 2.11 ). The arbitrary constant in ( [3.1| ) has been determined by imposing 
the condition that at z = we must have p = p M . Similarly, solving for the ^-dependent 
A as a function of the reshift, with the initial condition A(z = 0) = Ao, we find 

A(z; is)=A + P° M f(z) + p° c g(z) , (3.3) 

with 



/(*) 



1 



(l + z 



.3(1-!/) 



1 



9{z) 



z(z + 2) 



+ 



(1 + *) 



3(1-1/) 



1 



(3.4) 



(3.5) 



l-3v { 2 \-v 
Furthermore, since we have already obtained the functions p(z; v) and A = A(z; z/), we can 



use them in (L5) to get the explicit ^-dependent Hubble function H(z; v). It reads 



H 2 {z;v) = Hi 



+- 



[l + z 



,3(1-*/) 



1 



1 



- ft 



{l + zf 



l-2u 



(l + z 



,3(1-!/) 



(3.6) 
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For flat universes (|3.6|) shrinks to 



H 2 (z;u) = Hi 



(1 + *) 



3(1-*) 



M 



1 



(3.7) 



whereas the standard result in this case is I 



H\z) = Hi{l + n° M [(l + z) 8 -!]} 



(3. 



which is indeed recovered from Eq. (|3.7| ) in the limit v = 0. 

Notice that the kind of general behaviour that we have found for the CC evolution is 
of the form: 

A(t) = A + Z[H 2 (t)-H\t )]M 2 P , (3.9) 



where the coefficient £ is proportional to v: 



3^ 
8^r 



(3.10) 



Let us point out that many similar phenomenological equations for the CC evolution have 
been tried in the literature | 40| , 41], but the interesting thing in our particular realization is 



that it is motivated by the Renormalization Group, and therefore the structure (3J3) may 
emerge from a fundamental QFT formulation. 

Before applying our model to test the status of the present universe, it is important 
to check whether the model can be made compatible with nucleosynthesis. In fact, a non- 
vanishing v may have an impact not only in the matter-dominated (MD) era, but also in 
the radiation-dominated (RD) epoch. To this end we recall that in the RD era Eq. ( |1,7| ) 
must include the p ^ term. For photons the radiation density p R is related to pressure 
through p = (1/3) p R . With this modification we may solve again the differential equations 
in the radiation era. In this case it is more natural to express the above result in terms of 
the temperature. For the radiation density we find 



-4 I T o 



Pr(T] v) = —g*T ( — 



lis 



+ 



2-4i/ 



I) 



4(1-*) fT 



(3.11) 



where To ~ 2.75 K = 2.37 x 10 4 eV is the present CMB temperature. Of course for v — > 
we recover the standard result 



\ 4 



(3.12) 



with = 2 for photons and <?* = 3.36 if we take neutrinos into account. For the CC, 

A R (T; u) = A + p° M f R (T) + p° c g R (T) , (3.13) 

with 



f R (T) 



1 - v 



r \ 4(1-1/) 

5b) ^ 



(3.14) 
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and 



9R(T) 



k \T 2 -T, 



2-4i/l T 2 



r \ 4(1-1/) 

5b. 



(3.15) 



When comparing the relative size of the CC, Eq. ( [3.13 ), versus the radiation density, Eq 



(3.11), at the time of the nucleosynthesis we naturally require that the former is smaller 

npose 

~|i/|<1. (3.16) 



than the latter. Therefore, for small k, we impose 

A*(T) 



Pr{T) 



4. Numerical results 



After the nucleosynthesis restriction on the i/-parameter, the question is whether there is 
still some room for useful phenomenological considerations at the present matter epoch. 
Fortunately, the answer is yes. A convenient fiducial value for the cosmological index v is: 

= * 2.6 x 1(T 2 , (4.1) 



which corresponds to M = Mp in ( 2.11| ). We will use this value, and some multiples of it 



(with both signs) to perform the numerical analysis p0|| . The maximum value that we will 
tolerate for v is v = 0.1, which still respects the nucleosynthesis bound ( p. 16 ). Let us also 



circumscribe the numerical analysis to the flat case, k = 0. The evolution of the matter 
density and of the CC is shown in Fig.|l]a,b. These graphics illustrate Eq. ( |3.1| ) and (|Q|). 
As a result of allowing a non-vanishing ^-function for the CC (equivalently, v ^ 0) there 
is a simultaneous, correlated variation of the CC and of the matter density. 

Comparing with the standard model case v = (see the Fig.|a,b), we see that for a 
negative cosmological index v the matter density grows faster towards the past {z — ► oo) 
while for a positive value of v the growing is slower than the usual (1 + z) 3 . Looking 
towards the future (z — > —1), the distinction is not appreciable because for all v the 
matter density goes to zero. The opposite result is found for the CC, since then it is for 
positive v that A(z; v) grows in the past, whereas in the future it has a different behaviour, 
tending to different (finite) values in the cases v < and < v < 1, while it becomes 
— oo for v > 1 (not shown). 

In the phenomenologically most interesting case \v\ < 1 we always have a null density 
of matter and a finite (positive) CC in the long term future, while for the far past yields 
A = ±oo depending on the sign of v. In all these situations the matter density safely 
tends to +oo. One may worry whether having infinitely large CC and matter density in 
the past may pose a problem to structure formation. From Fig.[l|a,b it is clear that there 
should not be a problem at all since in our model the CC remains always smaller than the 
matter density in the far past, and in the radiation epoch z > 1000 we reach the safe limit 



(3.16). Actually the time where A(z; v) and pm(z; v) become similar is very recent. 

We may consider another relevant exponent describing how the universe evolves, the 
deceleration parameter This one is fully sensitive to the kind of high-z SNe la data 
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Figure 1: (a),(b) Future and past evolution of the matter density pm{z]v) and the cosmological 
constant A(z; v) for a flat universe (k — 0) and for different values of the fundamental parameter 
v of our model (vq is defined in Eq.(4.1)). In both cases = 0-3 an d = 0.7. The solid line 



represents the standard model case (A = const.), whereas the various kinds of dots/dashes represent 
different amounts of evolution of A J|(| . 

under consideration. By definition, 

a H 1 . , 1 dH 2 (z: v) 



where the ^-dependent Hubble parameter H 2 (z; v) is given in Eq. ( |3.6[) . It is interesting to 
look at the deviations of the deceleration parameter with respect to the standard model. 
It is well-known that there are already some data on Type la supernovae located very near 



the critical redshift z* where the universe changed from deceleration to acceleration [42]. 
But of course the precise location of z* depends on the FLRW model and variations thereof. 
In our case z* should depend on our cosmological index v, i.e. z* = z*(v). For simplicity 
in the presentation, let us consider the flat case. The transition point between accelerated 
and decelerated expansion is a function of v : the more negative is v, the more delayed is 
the transition (closer to our time)- see Fig. ||. If v > 0, the transition occurs earlier (i.e. 
at larger z). While in the standard case, and for a flat universe, the transition takes place 
at redshift 

/ qo 

z* = -l+V2-^f -0.67, (4.3) 

it would have occurred at z = 0.72 and z = 0.78 for v = vq and v = 2vq respectively, and 
at z = 0.63 and z = 0.59 for v = -u and v = -2u (Cf. Fig.@). For v = (-0.1, +0.1) the 
effect is quite large, namely the transition would be at z = (0.53, 0.91) and hence there is 
a correction of (—21%, +36%) with respect to the standard case. 
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Figure 2: Evolution of the deceleration parameter, q, for a flat universe (k — 0) with £l" f = 0.3 and 
= 0.7. We see that the transition from a decelerated (q > 0) to an accelerated (q < 0) universe 
takes place earlier in time (larger redshifts) for larger v > (see text for the exact values) p3| . 



It should be clear that our approach based on a variable CC departs from all kind 
of quintessence- like approaches, in which some slow-rolling scalar field x substitutes for 
the CC. In these models, the dark energy is tied to the dynamics of the self-conserved 



X field; i.e. in contrast to (1.7) there is no transfer between x-dark energy and ordinary 



forms of energy. The phenomenological equation of state is defined by p x = w x p x . The 



term —2 A on the r.h.s. of Eq. (1.6) must be replaced by p x + 3p x = (1 + 3w x )p x . In 
order to get accelerated expansion in an epoch characterized by p = and p — » in the 
future, one must require — w- < w x < —1/3, where usually W- > —1 in order to have a 
canonical kinetic term for x 5 - The particular case w x = — 1 corresponds to a quintessence 
field exactly mimicking the cosmological constant term. In fact, for a cosmological term A, 
whether constant or variable, the only possible equation of state is the one corresponding 
to u>a = —1, as it obvious from the definition of T^ u in Eq. (fh^). Although p x and p x are 
related to the energy-momentum tensor of x> t ne dynamics of this field is unknown because 
the quintessence models do not have an explanation for the value of the CC. Therefore, the 
barotropic index w x is not known from first principles. In particular, one cannot exclude 
it may have a redshift dependence, which can be parametrized with two parameters as 
follows: 

^ = w x = w + wi(l - a) = w + wit-— ■ (4-4) 
Px t + z 

Finding a non-vanishing value of W\ implies a redshift evolution of the equation of state 
for the x field ||. For completeness, consider the modification on the Hubble parameter 



'One cannot completely exclude "phantom matter-energy" (w_ < —1) and generalizations thereof J 
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Distribution 


Data 


Prior au M 




o- u 


SNAP 

(1 year) 


50 SNe < z < 0.2 
1800 SNe 0.2 < z < 1.2 
50 SNe 1.2 < z < 1.4 
15 SNe lA<z< 1.7 


None 


0.1 


±0.10 


SNAP 


as above 


0.03 


0.1 


±0.06 


SNAP 


3 years 


None 


0.1 


±0.06 


SNAP 


3 years 


0.03 


0.1 


±0.04 




2000 SNe 0.2 < z < 1.7 


0.03 


0.1 


±0.05 


Distr.2 


250 SNe < z < 1 
1750 SNe 1< z < 2 


0.03 


0.1 


±0.02 



Table 1: Determination of v with SNAP data and with other two distributions. In all cases we 
assume a flat universe. When a prior on Qm and its error is added we use Qm = 0.3 ± 0.03 |3(J. 



introduced by quintessence models. In this case, and in the flat case, it is easy to see that 
H 2 (z;w ,wi) = Hi {n° M (l + z) 3 



+ (i-n° M )(i + z ) 



3(l+ioo +w\) 



exp 



-3 wi 



l + z 



(4.5) 



This equation reduces to the standard one ( |3.8| ) for w x = —1 (wq = —l,wi = 0), as 
expected. Comparison between ( fl~5| ) and ( p.7| ) can be useful to identify the differences 
between RG models and quintessence models of the dark energy. 

Fitting cosmological models to high-z supernovae data is usually performed via the so- 
called magnitude-redshift relation H, [| . One starts from the notion of luminosity distance, 
di related to the received flux T and the absolute (intrinsic) luminosity C through the 
geometric definition [jjj: 

C (4-6) 



Then the logarithmic relation between flux and the (theoretical) apparent magnitude reads 
m th (z,H ,n° M ,n° A ) = M + 51og 10 [H d L (z,H ,n M ,n A )] , (4.7) 
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where M is the absolute magnitude (believed to be constant for all Type la supernovae, 
assuming they are real "standard candles"), and 

M = M - 5 log 10 H + 25 = M - 5 log 10 h + 42.38 . (4.8) 

The model dependence is encoded in the luminosity-distance function d^ = di(z, Hq, f^, ^ A ), 
given in our case by p(| 



iL( *MM:) - * fr/m jf , (4.0) 



with 



sinx, O^- < 

x, n° K = (4.10) 

sinhx, > • 



Here the expansion rate H(z,n° M ,n° A ;u) is given by Eq. (|3.6|) . In the flat case it reduces 
to the two-parameter function ( |3.7D , and in quintessence models to the three-parameter 
function fl4.5|) . Usually a prior on Q° M (from cluster dynamics) can be accepted (e.g. O^- = 
0.3), which narrows down the number of parameters to one (RG) and two (quintessence). 

We can use the magnitude-redshift relation defined above to test various distributions, 
including the one foreseen by SNAP . In order to determine the cosmological parameters 
we use a x 2 -statistic test, where \ 2 is defined by the difference between the theoretical 



apparent magnitude and the observed one (see Ref. [3C] for details). The idea is that the 
difference between models grows with redshift. The existing sample of SN la data is amply 
compatible with v ^ 0, but it does not pin down a narrow interval of values for this 
parameter. Trying, however, over several large (simulated) distributions of supernovae, 
and marginalizing over Ai , one obtains the results shown in Table ||. Distribution 1 is very 
similar to the SNAP one but with most of the data homogeneously distributed between 
z = 0.2 and z = 1.7. Distribution 2 extends data up to redshift z = 2. We see that we can 
determine v to within ±(20 — 60)% for v = 0.1, depending on the distribution. (Smaller 
values of v imply smaller precision.) The situation is similar to the determination of the 
evolution of the equation of state, Eq. ( |4.4[) . Indeed, if one performs a general (model- 
independent) fit of the present SN la data to quintessence models, leaving free the two 
parameters wo and w\ in Eq. (fO|), one finds that the values wq > —1/3 (decelerated 
universe) are ruled out at a high significance level (for Q^j < 0.4), thereby supporting the 



existence of dark energy. Nevertheless, the very same fit is highly insensitive to w\ [39 
[l(|. On the other hand SNAP will be able to determine {£1® m ,wq) to within small errors 
(3%, 5%), and will significantly improve the determination of the time- variation parameter 
u>i, but only up to 30% at most (each parameter being marginalized over others) [p5| . 
Moreover, the degeneracies among their combinations weaken substantially the tightness 
of these bounds. 
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5. Concluding remarks 

We have considered a Friedmann-Lemaitre-Robertson- Walker model with time-evolving 
cosmological term: A = A(i). The evolution of the CC is due to quantum effects that 
can be described within the Renormalization Group (RG) approach. The RG scale /x is 
identified with the Hubble parameter H at the corresponding epoch, as first proposed 



in [21]. Although the ^-function for A is just proportional to the fourth power of the 
masses in the UV regime, the phenomenon of decoupling in QFT leads to an inverse power 
suppression by the heavy masses at low energies (IR regime). Thus, in the present day 
universe, one may expect a modified RG equation characterized by a "soft decoupling" 
behaviour j3\ ~ H 2 M 2 JH], ||, ||, |o|. The effective mass scale M summarizes the 



presence of the heavy degrees of freedom. This peculiar form of decoupling can be envisaged 
from: i) the Appelquist-Carazzone (AC)-decoupling theorem, ii) general covariance of the 
effective action, and also from iii) the non-fine-tuning hypothesis on the n = 1 terms of 
Pa ~ Cf. Eq. (|2.5[)-(fO|)- which insures that the coefficient of the quadratic contribution 
H 2 M 2 does not vanish. This particular form of decoupling is a specific feature of the CC 
because it is of dimension four. There is no other parameter either in the SM or in GUT 
models with such property. 

In constructing this semiclassical "RG-cosmology" we have explored the possibility 
that the heaviest d.o.f. may be associated to particles having the masses just below the 
Planck scale, hence M is of order Mp. This assumption is essential to implement the soft 
decoupling hypothesis within the \i = H setting, for the present value of HqM 2 is just of 
the order of the CC. This fact insures a smooth running of the cosmological term around 
the present time. In this model the /3A-function has only one arbitrary parameter v ( 2.11| ) 



proportional to the dimensionless ratio M 2 /Mp, and as a result the model has an essential 
predictive power. In general we expect \v\ <C 1 from phenomenological considerations, 
mainly based on the most conservative hypotheses on nucleosynthesis. 

As the variation of the CC is attributed, in this model, to the "relic" quantum effects 
associated to the decoupling of the heaviest degrees of freedom below the Planck scale, a 
time dependence of the CC may be achieved without resorting to scalar fields mimicking 
the cosmological term ( "quintessence" ) or to modifications of the structure of the SM of the 
strong and electroweak interactions and/or of the gravitational interactions. This proposal, 
therefore, offers an excellent opportunity to explore the existence of sub-Planck physics in 
direct cosmological experiments, such as SNAP (and the very high-z SNe la data to be 
obtained with HST). For v < 0.1 corrections to some FLRW cosmological parameters 
become as large as 50% or more, which could not be missed by these experiments. 

Whether this RG-cosmology can be easily distinguished from quintessence models re- 
quires further considerations along the lines of previous studies on evolving dark energy 
p9| , [L6|]. The present model, however, can elude some of the difficulties (related to the 
degeneracies of the kinematical and geometrical measurements fifl), in that we predict 
not only an evolving vacuum energy, but also a correlated (^-dependent) departure of the 
matter-radiation density from the standard model prediction. From a more fundamental 
point of view, the sole fact that our FLRW scenario with running cosmological constant 
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can compete on the same footing with quintessence models shows that standard quantum 
field theory in curved space-time may contain the necessary ingredients from which one 
can build up a time-evolving cosmological term without need of artificial ("just so") scalar 
fields. Both types of models can be thoroughly checked by the SNAP and upgraded HST 
experiments If these experiments will detect the redshift dependence of the CC similar 
to that which is predicted in our work, we may suspect that some relevant physics is going 
on just below the Planck scale. If, on the contrary, they unravel a static CC, this may im- 
ply the existence of a desert in the particle spectrum below the Planck scale, which would 
be no less noticeable. In this respect let us not forget that the popular notion of a GUT 
(perhaps in the form of string physics) near the Planck scale remains, at the moment, as 
a pure (though very much interesting!) theoretical speculation, which unfortunately is not 
supported by a single piece of experimental evidence up to now. Our framework may al- 
low to explore hints of these theories directly from astrophysical/cosmological experiments 
which are just round the corner. If the results are positive, it would suggest a direct link 
between the largest scales in cosmology and the shortest distances in high energy physics. 
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